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[l]  Ocean  predictability  skill  is  investigated  using  a  Gulf 
of  Mexico  nowcast/forecast  model.  Power  law  scaling  is 
found  in  the  mean  square  error  of  displacement  between 
drifting  buoy  and  model  trajectories  (both  at  50  m  depth). 
The  probability  density  function  of  the  model  valid 
prediction  period  (VPP)  is  asymmetric  with  a  long  and 
broad  tail  on  the  higher  value  side,  which  suggests  long¬ 
term  predictability.  The  calculations  demonstrate  that  the 
long-tenn  (extreme  long  such  as  50-60  day)  predictability 
is  not  an  "outlier"  and  shares  the  same  statistical  properties 
as  the  short-term  predictions.  INDEX  TERMS:  4263 
Oceanography:  General:  Ocean  prediction;  3210  Mathematical 
Geophysics:  Modeling;  4532  Oceanography:  Physical:  General 
circulation;  4512  Oceanography:  Physical:  Currents 

1.  Introduction 

[2]  Two  kinds  of  predictability  exist  in  ocean  (and 
atmospheric)  models.  The  predictability  of  the  first  kind 
(second)  describes  the  model  performance  due  to  imprecise 
initial  conditions  (boundary  conditions  or  model  parame¬ 
ters)  [Lorenz,  1984],  For  a  perfect  model,  a  universal 
exponential  law  is  found  in  the  first  type  predictability  for 
prediction  error  growth  during  the  initial  stage  (linear 
regime)  due  to  infinitesimal  and  small  initial  errors  [Lorenz, 
1984;  Robinson,  1996]. 

[3]  Neither  the  model  dynamics  nor  the  model  forcing  is 
identical  to  that  of  the  real  ocean.  For  example,  the  lateral 
transport  through  straits  and  fluxes  between  the  atmosphere 
and  ocean  are  not  precisely  known  apriori.  Uncertain  model 
parameters  and  unresolved  spatial  and  temporal  scales 
distort  the  dynamical  equations  (the  evolution  law).  The 
forecast  sensitivity  to  such  distortions  is  referred  to  as  the 
predictability  of  the  second  kind  [Lorenz,  1984],  Chu  [1999] 
has  investigated  the  two  kinds  of  predictability  using  the 
Lorenz  attractor. 

[4]  On  the  base  of  the  first-passage  time  (i.e.,  the  time 
period  for  a  moving  particle  inside  a  domain  to  first 
impacts  its  boundary),  Chu  et  al.  [2002a,  2002b]  proposed 
that  the  valid  prediction  period  (VPP)  effectively  repre¬ 
sents  predictability.  VPP  is  defined  as  the  time  period 
when  the  prediction  error  first  exceeds  a  pre-determined 
criterion  e  (i.e.,  the  tolerance  level).  The  conditional 
probability  density  function  of  VPP  with  a  given  initial 
error  satisfies  the  backward  Fokker-Planck  equation 
(called  Pontryagin-Kolomogorov  equation  in  Russian  liter¬ 
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ature).  Using  VPP  as  a  quantitative  measure  for  prediction 
skill,  both  linear  and  nonlinear  regimes  of  forecast  errors 
were  found  in  the  low-order  atmospheric  Lorenz  model 
[Lorenz,  1984].  In  the  linear  regime  (usually  for  the 
predictability  of  the  first  kind  during  the  initial  stage), 
the  error  growth  follows  the  exponential  law  (linear 
scaling  law)  with  VPP  represented  by  the  e-folding  time. 
In  the  nonlinear  regime,  the  error  growth  does  not  follow 
the  exponential  law,  and  VPP  cannot  be  represented  by  the 
e-folding  time. 

[5]  The  existence  of  the  power  law  (nonlinear  scaling 
law)  is  well  known  in  oceanic  and  atmospheric  transport 
processes  [Osborne  et  al.,  1989;  Meyers  et  al.,  1994; 
Kyong-Hwan  and  Bowman,  2000],  Questions  arise:  Does 
the  power  law  exist  in  the  error  growth?  Where  does  it 
occur?  Since  the  exponential  law  was  found  in  the  predict¬ 
ability  of  the  first  kind,  it  is  reasonable  to  check  if  the  power 
law  of  error  growth  exists. 

[6]  The  goal  of  this  study  is  to  demonstrate  the  existence 
of  the  power  law  of  error  growth  using  a  Gulf  of  Mexico 
nowcast/forecast  model  and  the  drifter  buoy  data  and  to  find 
a  long-term  (50-60  day)  predictability.  The  predictability 
skill  is  found  to  have  a  long  tail  (on  the  high  value  side) 
following  the  power  decay  law  (t_x)  for  long  timescales, 
where  A  is  the  scaling  exponent.  It  is  well-known  [Karney, 
1983]  that  such  a  behavior  leads  to  a  slower  error  increase 
than  the  exponential  growth.  Thus,  the  error  growth  follow¬ 
ing  the  power  law  can  be  referred  to  as  the  long-term 
predictability. 

2.  Quantitative  Measures  of  Model 
Predictability  Skill 

[7]  Let  y(t)  =  [ydxj,  t),  y2(x2,  t),  . . .,  yN(xN,  t)]  be  a  set  of 
trajectories  from  N  drifter  buoys,  and  y (m)(t)  =  [j^m)(Pi ,  t), 
y2m)(p2,  t), . . .,  Jn'^Pn,  0]  t>e  a  set  of  model  trajectories  from 
N  synthetic  particles.  Here  (X| ,  x2 ,. .  .,xN)  and  (pi,  p2,  . . ., 
pN)  denote  the  initial  positions  of  drifter  buoys  and  syn¬ 
thetic  particles,  respectively.  Thus,  the  forecast  error  is 
represented  by 

z(r)  =  yW(t)-y(t), 

z(f0)  =  Z0  =  [Pj  -  Xl,p2  -  X2,  .  .  .  ,Pjv  -  Xjv]- 
The  error  mean  and  variance  are  computed  by 

L\  =  (z),i2  =  ((z  -  (z))'(z  —  (z)))  (1) 

where  the  bracket  indicates  the  ensemble  averaging.  If  the 
error  growth  follows  an  exponential  law 

L\  oc  eat,  L2  oc  ew<,  (2) 
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a  and  oj  are  the  growth  rate,  and  if  the  error  growth  is 
described  by  a  power  law 

L\  oc  ta,  L2  oc  ft,  (3) 

a  and  (3  are  scaling  exponents. 

[s]  The  power  law  (3)  can  be  interpreted  in  terms  of 
anomalous  transport  and  diffusion  [Weeks  etal. ,  1996].  If  the 
mean  displacement  grows  linearly  with  time  (a  =  1),  there  is 
a  well-defined  transport  velocity.  The  condition  a  7^  1  is 
referred  to  as  anomalous  advection,  and  the  condition  (3^1 
is  referred  to  as  anomalous  diffusion  with  0  <  (3  <  1  as  sub¬ 
diffusion,  1  <  (3  <  2  as  supper-diffusion,  and  (3  =  2  as  ballistic 
dispersion  occurring  in  the  shear  flows. 

[9]  VPP  is  represented  by  a  time  period  t  at  which  |z| 
first  exceeds  a  predetermined  tolerance  level  e.  The  prob¬ 
ability  density  function  (PDF)  of  VPP,  W(l,  z0,  e),  with  the 
initial  error  z0  and  tolerance  level  6,  satisfies  the  backward 
Fokker-Planck  equation  [Chu  el  al.,  2002a,  2002b;  Chu 
and  Ivanov,  2002],  The  predictability  skill  can  also  be 
verified  using  the  probability  density  function  of  prediction 
(PDP),  which  is  defined  by  successful  prediction  for  the 
period,  t  —  t0  [Ivanov  et  al.,  1994;  Ivanov  and  Margolina, 
1999], 

PDP(t  >  t-  t0,z0,e)  =  P(t0,  z0,e,t-  t0) 

‘-to 

=  1  ~  J  W(r,z0,£)dT,  (4) 

0 

where  to  is  the  initial  time.  If  a  power  law  (3)  exists  for  the 
prediction  errors,  PDP  should  also  follow  the  power  law, 

P(to,  Zo,  e,  t  -  to)  ~  r7  for  large  t.  (5) 


3.  The  Gulf  of  Mexico  Nowcast/Forecast 
System  and  Drifter  Data 

[10]  The  Gulf  of  Mexico  real  time  nowcast/ forecast 
system  is  taken  as  an  example  to  show  the  existence  of 
the  power  decay  law  and  the  range  of  the  three  power 
parameters  (a,  (3,  y,)  in  the  predictability  skill.  This  now¬ 
cast/forecast  system  [Kantha  et  al,  1999;  Schaudt  et  al., 
2001;  Kantha  and  Clayson,  2000]  is  based  on  the  Univer¬ 
sity  of  Colorado  version  of  the  Princeton  Ocean  Model 
(POM)  with  1/12°  horizontal  resolution.  Real-time  altimet- 
ric  sea  surface  height  (SSH)  anomalies  derived  from  NASA/ 
CNES  TOPEX  and  ESA  ERS-1/2  altimeters,  and  composite 
MCSST  data  derived  from  NOAA  AVHRR  assimilated  into 
the  model  in  a  continuous  data  assimilation  mode  are  used 
to  produce  a  nowcast  and  a  four-week  forecast.  Kantha  and 
Clayson  [2000]  found  that  forecast  retains  considerable  skill 
to  about  1-2  weeks,  beyond  which  the  forecast  begins  to 
deviate  from  reality. 

[11]  The  Gulf  of  Mexico  velocity  data  at  50  m  depth 
were  archived  every  six  hours  from  the  nowcast/forecast 
system  for  six  months  beginning  on  July  9,  1998. 
The  observational  data  were  obtained  for  50  satellite- 
tracked  drift  buoys  drouged  at  50  m  depth.  These  buoys 
were  manufactured  by  the  Applied  Technology  Associates, 
deployed  by  the  Florizon  Marine  Inc.,  and  tracked  by 


Figure  1.  Statistical  characteristics  of  VPP  for  zero  initial 
error  and  55  km  tolerance  level:  (a)  histogram  of  VPP,  (b) 
PDF  of  VPP  computed  using  the  Epanechnikov  kernel 
density,  (c)  PDP  for  VPP  between  0  to  20  days,  and  (d)  PDP 
for  VPP  between  5  to  17  days. 

the  Louisiana  State  University’s  Coastal  Studies  Institute. 
Navigation  errors  for  the  buoy  position  are  less  than 
100  m. 

4.  Power  Law  of  Error  Growth  in  Non- 
Ensemble  Prediction 

[12]  Since  the  number  of  drift  buoys  is  limited  (50  in  this 
study),  the  point  stochastic  process  is  used  to  create  more 
realizations  to  calculate  the  ensemble  mean  of  the  predict¬ 
ability  skill  [Tikhnov  and  Khimenko,  1998].  Moving  along 
the  drifter  trajectory  y(t),  points  are  randomly  selected  and 
the  predictability  skill  is  calculated  for  each  of  these  points. 
Such  a  process  is  called  the  point  stochastic  process.  An 
average  over  these  points  is  equivalent  to  the  ensemble 
average  due  to  the  ergodic  feature  of  the  trajectories 
[Dymnikov  and  Filatov,  1977;  Tikhonov  and  Khimenko, 
1998], 

[13]  VPP  is  computed  for  non-ensemble  prediction  with 
tolerance  level  (e)  of  55  km  and  without  initial  error 
(z0  =  0).  Its  PDF,  W(t,  0,  55  km)  (Figure  lb),  calculated 
using  the  Epanechnikov  kernel  density  [Good,  1996]  from 
the  histogram  (Figure  la),  shows  a  non-Gaussian  distri¬ 
bution  with  a  narrow  peak  and  a  long  tail  in  the  domain 
of  long-term  prediction.  Successful  10-15  day  predictions 
(Figure  lc)  are  abnormally  long  (called  extreme  long 
predictions)  compared  to  the  mean  VPP  (3.2  days)  and 
the  most  probable  VPP  (2.4  day).  The  tail  of  PDP 
follows  the  power  law  (Equation  5)  with  the  power 
exponent,  y  =  2.2  (Figure  Id). 

[14]  The  error  mean  ( L\ )  and  variance  (L2)  are  calcu¬ 
lated  at  1320  points  generated  using  the  stochastic  point 
process.  The  error  growth  from  July  9,  1998  (day-1)  to 
August  3,  1998  (day-25)  shows  the  existence  of  the 
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predictability  of  the  first  kind  (following  an  exponential 
law). 

6.  Long  Term  Predictability  With  Various 
Tolerance  Levels 

[16]  To  understand  the  effect  of  tolerance  level  on  the 
long-term  predictability  (or  the  power  law),  the  PDF  of  VPP 
is  computed  with  zero  initial  error  and  with  varying  e 
(Figure  4).  All  the  PDFs  show  a  similar  pattern  to  PDF 
for  £  =  55  km  (Figure  lb):  a  non-Gaussian  distribution  with 
a  narrow  peak  and  a  long  tail  in  the  domain  of  long-term 
prediction.  The  larger  the  tolerance  level  (e),  the  longer  the 
tail.  Variation  of  tolerance  level  does  not  distort  the  power 
law  and  the  long-term  prediction.  Extreme  long-term  pre¬ 
diction  exists  even  for  short  mean  and  most  probable  VPP. 
For  example,  the  VPP  of  extreme  long-term  prediction  can 
reach  50-60  days  for  8  >  165  km  when  the  mean  VPP  is 
within  10-15  days  (Figure  4b).  Moreover,  the  extreme 
long-term  predictions  are  not  "outliers"  [Lvov  et  al., 
2001],  since  the  same  PDF  is  applied  to  short-  and  long¬ 
term  predictions. 

7.  Statistical  Similarity 

[17]  The  physical  reason  for  the  existence  of  extreme 
long-term  prediction  is  the  statistical  similarity  between 
modeled  and  observed  circulation  patterns.  To  understand 


Figure  2.  Scaling  behavior  of  the  mean  error  ( L , )  growth 
for  initial  error  level  of  (a)  zero,  (b)  2.2  km,  and  (c)  22  km, 
respectively. 


power  laws  in  L\  and  L2  with  the  scaling  exponents  a  = 
0.72  ±0.02  (Figure  2a)  and  [3  =  1.34  ±0.02  (Figure  3a), 
respectively. 

5.  Power  Law  of  Error  Growth  in  Ensemble 
Prediction 

[15]  For  ensemble  prediction,  the  initial  position  of  syn¬ 
thetic  particles  is  perturbed  from  the  initial  position  of 
drifter  buoys 


pk  =  xk  +  rk,  (6) 

where  rk  is  the  initial  perturbation  and  its  magnitude 
varies  from  2.2  to  22  km.  The  error  mean  (Lx)  and 
variance  ( L2 ),  obtained  from  the  ensemble  prediction,  also 
follow  the  power  law.  For  example,  the  power  for  Lt 
varies  from  a  =  0.71  ±  0.02  for  initial  error  level  of  2.2 
km  (Figure  2b),  to  a  =  0.63  ±  0.03  for  initial  error  level  of 
22.5  km  (Figure  2c).  The  power  for  L2  varies  from  (3  =  1.4 
±  0.02  for  initial  error  level  of  2.2  km  (Figure  3b),  to  [3  = 
1.3  ±  0.02  for  initial  error  level  of  22.5  km  (Figure  3c). 
Ensemble  prediction  (uncertainty  inserted  in  initial  condi¬ 
tion)  does  not  distort  the  power  law.  It  only  changes  the 
value  of  scaling  exponents  a  and  (3.  Thus,  the  forecast 
errors  for  the  predictability  of  the  second  kind  (following  a 
power  law)  grow  considerably  slower  than  for  the 


Figure  3.  Scaling  behavior  of  the  error  variance  (L2) 
growth  for  initial  error  level  of  (a)  zero,  (b)  2.2  km,  and  (c) 
22  km,  respectively. 
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Figure  4.  PDF  of  VPP  calculated  with  different  tolerance 
levels:  (a)  from 27.5  to  110  km,  and(b)  from  137.5  to  2 10  km. 

such  a  similarity,  empirical  orthogonal  function  (EOF) 
analysis  is  conducted  for  the  Gulf  of  Mexico  circulation. 
The  EOFs  {''I',}  and  the  corresponding  amplitudes  [/!,(/)! 
for  the  streamfunction  are  computed  from  the  model  veloc¬ 
ity  fields  at  50  m  depth  during  the  six  month  period  (April 
10-October  10,  1998)  using  the  method  proposed  by 
Penenko  and  Protasov  [1978],  The  first  ten  EOFs  consist 
of  more  than  90%  of  the  total  variability. 

[is]  The  velocity  along  the  drift  trajectories  is  represented 
by 

10 

V  =  ^5,(t)k-  V'T;  (7) 


where  k  is  the  unit  vector  in  the  vertical  (upward  positive), 
and  V  is  the  horizontal  gradient  operator.  The  mode 
amplitudes  {Bi(t), . .  ,,Bl0(t)}  are  estimated  from  the  drifter 
buoy  data  [Eremeev  et  al.,  1992],  The  sporadic  velocity 
vector  field  of  the  drifter  buoys  at  00:00  on  July  9,  1998 
(Figure  5)  is  reconstructed  into  a  gridded  velocity  vector 
field  (Figure  5b)  using  Equation  (7).  The  modeled  velocity 
vector  field  (Figure  6a)  is  consistent  with  the  reconstructed 
velocity  vector  field  (Figure  6b).  It  is  noticed  that  the 
domain  in  Figures  5  and  6  is  at  50  m  depth. 

[19]  The  two  curves  {/t,(t)}  and  {S,(t)}  in  any  amplitude 
pair  (i  =  1,  2,. . .,  10)  have  seven  crossing  points  within  the 
6  month  period.  At  the  crossing  points,  there  is  no  differ¬ 
ence  between  the  modeled  and  reconstructed  (from  drift 
buoy  data)  velocity  fields.  The  Kolmogorov-Smirnov 
probabilities  for  {A  1,  5.},  \A2,  Bo}, . . .,  {A10,  Bl0}  are, 
0.63  ±  0.4,  0.79  ±  0.3,. . .,  0.39  ±  0.2  respectively.  This 
confirms  the  statistical  similarity  between  {Aj( t)}  and 

m)i 


Figure  5.  Velocity  vectors  of  drifters  calculated  on  00:00 
July  9,  1998  at  50  m  depth. 


18 


86  “W  W“W  88  "W  84'W 

Figure  6.  Velocity  vectors  on  00:00  July  9,  1998  obtained 
using  (a)  the  Gulf  of  Mexico  model  and  (b)  the 
reconstruction  scheme  on  the  drift  data. 
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8.  Conclusions 

[20]  A  power  law  exists  in  ocean  model  error  growth.  The 
probability  density  function  of  VPP  for  the  Gulf  of  Mexico 
nowcast/forecast  system  is  asymmetric  with  a  long  and 
broad  tail  on  the  higher  value  side,  which  indicates  long¬ 
term  (50-60  day)  predictability.  Such  a  long  tail  corre¬ 
sponds  to  a  power  law  for  the  error  growth.  Individual 
model  predictions  can  be  valid  for  a  long  period  (long-term 
prediction),  even  for  short  mean  and  the  most  probable  VPP 
(10-15  day).  The  extreme  long  and  short  predictions  share 
the  same  statistics. 
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